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, . . Abstract 

X 

ci I We discuss the Sakai-Sugimoto model at finite temperature and finite 

chemical potential. It is a holographic model of large Nc QCD with Nf 
massless quarks based on a D4/D8-D8 brane system. The near horizon 
limit of the D4-branes and the probe approximation of the D8-D8 pairs 
allow us to treat the D4-branes as a gravitational background and the D8- 
D8 pairs as a probe which does not affect the background. We propose 
that the asymptotic value of a U(l) gauge field on the D8-D8-branes is 
identified with the chemical potential for the baryon number. Using this 
chemical potential we analyze the phase structure of this model and find 
a chiral symmetry phase transition of the first order. 
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1. Introduction 



The AdS/CFT correspondence [H [21 [3] (see [1] for a review) is a useful duality 
between a string theory in [d + l)-dimensional anti de Sitter spacetime (times a 
compact space) and a d-dimensional conformal field theory. The AdS/CFT corre- 
spondence can be extended to more general cases of the string/gauge duality for 
non-conformal and non-supersymmetric theories. In this scheme one can discuss 
some features of the low energy QCD such as confinement and spontaneous chiral 
symmetry breaking. Such an approach to low energy behaviors of QCD in terms of 
the string/gauge duality is often called the holographic QCD |5l [6], O [H O [10] (and 
references therein). 

One of the interesting recent developments in the holographic approach to QCD 
is the D4/D8-D8 model proposed by Sakai and Sugimoto [HI [12]. This brane system 
consists of Nc D4-branes compactified on S^ and probe Nf D8-D8-brane pairs. The 
D4-branes are described by the extremal brane solution in the near horizon limit 
with one of the spatial directions along its world- volume compactified on S^. This 
gravitational background is dual to a five-dimensional gauge theory, which looks 
four-dimensional at energy scale below the compactification scale. This description 
is valid for the case 1 <C Qym^c "C l/gyM^ where gvM is the four-dimensional 
gauge coupling. Imposing periodic boundary conditions on the bosons and anti- 
periodic ones on the fermions along the compactified direction, supersymmetry is 
explicitly broken. The scalars and the fermions on the D4-branes become massive 
and are decoupled from the system at low energy. Thus one obtains a U(A'"c) pure 
gauge theory. To describe quarks in the fundamental representation of the gauge 
group U{Nc) one introduces Nf D8-D8 pairs into the D4 background. The probe 
approximation Nf Nc ^ allows us to treat the Nf D8-D8 pairs as a probe, 
which does not affect the D4 background. A string connecting the D4-branes and 
the D8-branes (D8-branes) represents a massless left-(right-)handed quark with Nf 
flavors. Therefore one obtains a four- dimensional U{Nc) gauge theory with Nf 
flavored massless quarks in the fundamental representation of the gauge group at 
low energy. The U(A^/)l x U(A^/)/j symmetry of the D8 and D8-branes represents a 
chiral symmetry of the quarks. It was shown in ref. [11] that this chiral symmetry 
is broken to a diagonal symmetry U(A^/)y because of a configuration of the probe 
D8-D8 pairs. By the topology of the D4-brane geometry they must have a curved 
configuration in which the D8-branes and the D8-branes are connected each other 
(see Fig. [2]), which breaks the chiral symmetry. 
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The Sakai-Sugimoto model was also studied at finite temperature T [T3l fTH ITS]. 
There are two D4-brane geometries which represent a low temperature phase and 
a high temperature phase respectively. A phase transition between these phases 
occurs at a certain critical temperature T^. This transition corresponds to a con- 
finement /deconfinement transition in a dual gauge theory [16]. For each of the 
phases one can introduce probe D8-D8 pairs. The only configuration of the D8-D8 
pairs which can realize in the low temperature phase T < Tc is a curved D8-D8 
configuration as in the zero temperature case. Thus the chiral symmetry is always 
broken in the low temperature phase. On the other hand, another configuration 
exists in addition to the curved one in the high temperature phase T > Tc. The 
new configuration consists of straight D8-branes and D8-branes disconnected each 
other (see Fig. [3]). The chiral symmetry is unbroken for this configuration. A chiral 
symmetry phase transition can occur between these two configurations at a certain 
temperature [131 [H]. 

The purpose of the present paper is to analyze the Sakai-Sugimoto model at 
finite temperature T and finite baryon number chemical potential /i. We introduce a 
non-vanishing background U(l) gauge field on the probe D8-D8-brane world- volume. 
The asymptotic value of this gauge field background is related to the baryon number 
chemical potential. There are several related works in which chemical potentials are 
introduced as asymptotic values of gauge fields [TTJ dSl [ISl EOl EI]- The gauge 
fields and the chemical potentials considered there, however, are those for the R 
symmetry in the bulk geometry or for the isospin symmetry on the probe brane 
world- volume. In contrast, we consider the chemical potential for the baryon number 
U(l)v symmetry on the probe brane. The D8-D8-brane configuration and the gauge 
field background are determined by equations derived from the Dirac-Born-Infeld 
effective action on the world-volume of the probe branes. By solving these equations 
and comparing values of the effective action for the solutions we discuss a chiral 
symmetry breaking as in refs. [TH [131 E] ■ 

We can summarize our results as follows. In the low temperature phase T < Tc 
of the confinement / deconfinement transition there is a unique solution for the probe 
brane configuration and the gauge field background. As in the case without the 
gauge field background the probe branes have a curved configuration and the chiral 
symmetry is always broken. In the high temperature phase T > there are two 
types of solutions. One solution has a curved brane configuration and the chiral 
symmetry is broken. The other solution has a straight brane configuration and the 
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Figure 1: The phase diagram of the dual gauge theory. 

chiral symmetry is unbroken. A phase transition between these two solutions occurs 
at certain temperature and chemical potential. The transition is of the first order. 
There is a critical value of the chemical potential above which the phase transition 
never occurs for any temperature. A qualitative feature of the phase diagram is 
shown in Fig. [H 

The organization of this paper is as follows. In section 2 we review the bulk 
geometry of the Sakai-Sugimoto model at finite temperature and set up the effective 
action for the probe brane configuration and the U(l) gauge field. In sections 3 and 
4 we obtain solutions for the brane configuration and the gauge filed and discuss 
the chiral symmetry breaking in the low and high temperature phases respectively. 
Section 5 is devoted to some discussions. 

While preparing the manuscript of this paper, we have received a paper [22], in 
which a gauge field background on the probe branes is used to represent the baryon 
number chemical potential in studying hadronic matters in the Sakai-Sugimoto 
model. After submitting the manuscript of this paper for publication, we have 
received a paper [23] discussing the chiral phase transition in the D4/D8-D8 model, 
in which an error in the first version of our paper was pointed out. 



2. D4/D8-D8 brane system 

The Sakai-Sugimoto model [HI [12] is based on a D4/D8-D8 brane system con- 



4 



sisting of compactified Nc D4-branes and Nf D8-D8-brane pairs transverse to the 
S^. The brane configuration of the system is 

t x"^ T u e"^ 

D4 O O O O O - - - ^ ^ (2.1) 
D8-D8 oo o o — ooooo 

with r and ^^'s being coordinates of and respectively. The period of r is denoted 
as 5t = 271 /MxK- In the large Nc hmit and the near horizon hmit the D4-branes 
are described by a bulk background geometry, which is a classical solution of the 
type llA supergravity in ten dimensions. Assuming Nf -C Nc the D8-D8 pairs are 
treated as a probe which does not affect the bulk background. 

The finite temperature behavior of the Sakai-Sugimoto model was discussed in 
fi3[ HH |15]. The bulk background geometry is represented by a metric with a 
periodic Euclidean time coordinate tE = it ^ + StE in addition to the periodic 
r. The period of tE is the inverse temperature 6tE = 1/T. There are two such 
Euclidean solutions which have an appropriate asymptotic boundary behavior. One 
of them is the Euclidean version of the extremal D4-brane geometry compactified 
on with the metric 

ds' = (^^^ ' {dtl + 6,,dx'dx^ + f{U)dr') + ' (^-^ + U^dni^ , 

/(f/) = l-^, Ukk = Ir'MI^, (2.2) 

where dVL\ is the metric of and B? = TCQgNclg with Qs and Ig being the string 
coupling and the string length. The parameter Ukk must be related to M^k as 
above to avoid a singularity of the metric at f/ = Ukk- With this relation the t-U 
submanifold has a cigar-like form with a tip at f/ = Ukk- The dilaton and the 
RR 3-form C3 are given by 

= 9s (^) F, = dCs = ^64, (2.3) 

where £4 and V4 are the volume form and the volume of S^. The other solution is 
the Euclidean version of the non-extremal D4-brane geometry 

ds' = (^^) ' {fiU)dtl + 6,,dx^dx= + dr') + (^^^ ' (^-^ + U'dni^ , 

m) = 1 - ^, Ut= '-^R'T^ (2.4) 
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with the dilaton and the RR 3-form given in eq. fl2.3p . The parameter Ut must be 
related to T as above to avoid a singularity of the metric at U = Ut- The tE-U 
submanifold has a cigar-like form with a tip at f/ = Ut- It is obvious that these two 
backgrounds are related by interchanging r, Ukk and tE, Ut- 

It was shown [161 El [13] that the background (I2.2p is dominant at low tempera- 
ture, while (I2.4p is dominant at high temperature by comparing values of the Eu- 
clidean supergravity action for these backgrounds. A phase transition between these 
backgrounds occurs at the temperature for Ut = Ukk, i-e. Tc = Mkk/{'^tt)- This 
phase transition is of the first order and represents a confinement /deconfinement 
transition [16] . 

In refs. [T3l [Ti] Nf D8-D8 pairs were introduced as a probe in the backgrounds 
fl2.2l) . fl2.4p . The effective action of the D8-branes consists of the Dirac-Born-Infeld 
action and the Chern-Simons term 

Sns = TsJ (fx Tr ^det ( (^Miv + 27ra'FM^) - J TtF\ (2.5) 

where qmn and Fmn = Qm^n — On Am — i [Am, An] (M, = 0, 1, ■ ■ ■ , 8) are the 
induced metric and the field strength of the U(A^/) gauge field Am on the D8-branes. 
Tg is the tension of the D8-brane and a' = 11 is the Regge slope parameter. The 
effective action for the D8-branes has a similar form. The total effective action has 
a gauge symmetry 

\J{Nf)L X \]{Nf)R = ^\]{Nf)L X ^\]{Nf)R X \]{l)v X U(l)^, (2.6) 

where U(A^/)l and U(A^/)i? are symmetries of Nj D8 and D8-branes respectively. 
It was argued in ref. [TT] that this gauge symmetry corresponds to a flavor chiral 
symmetry of quarks. The total effective action can be written in the form (12. 5p with 
the integrations being over the whole of the D8-D8 world-volume. We use this form 
of the effective action in the following. 

In refs. [HI [13] the gauge fields on the probe branes are treated as fluctuations 
representing the hadron spectrum. In this paper we consider a background gauge 
field. We assume that only the Euclidean time component of the U(l) gauge field has 
a non-vanishing background. We will see that it corresponds to an introduction of 
the baryon number chemical potential. We use a physical gauge for D8-brane world- 
volume reparametrizations and use the spacetime coordinates other than r as the 
world-volume coordinates. Then, D8 and D8-brane configurations are determined 
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by r as a function of those world-volume coordinates. We make an ansatz that Aq 
and r depend only on the coordinate U 

T = r{U), Ao = Ao{U). (2.7) 

By this ansatz the Chern-Simons term in eq. (12. 5p vanishes and does not concern 
us. 



3. Low temperature phase 



In the low temperature phase the geometry ( 12. 2p is dominant. Using the ansatz 
(12.71) the induced metric qmn on the probe D8-branes is given by 



+ 



du'+{^] u'dnl (3.1) 



where r' — '^'^ 



^jj. Then, the effective action of the D8-branes (12. 5p becomes 

3 



g N/TsYl I d'xdUU' 



9s 



/(r')^+^ (r^-(2WA^)^) 



(3.2) 



where A'q = ^. 



This action leads to equations of motion for t{U) and Ao{U) 



d 

dU 




0, 



'firf + {§Y{f-^-i2na'A',f) 
which can be easily integrated once. We obtain 

(^ {U)) 



(3.3) 



{2na'A',{U)y 



f{Uf [{^f (U^fiU) - U^fiUo)) + [f{U) - f{Uo) 



{§f (U^fiU) - UlfiU,)) + (f{U) - f{U,) (t) 



(3.4) 



Figure 2: A D8-D8-brane configuration in the low temperature phase. 

where C and Uq are integration constants. As in the zero temperature case [TT] and 
the low temperature phase in ref. [I3] we have imposed a condition t'{Uq) = oo. A. 
typical configuration of t{U) is shown in Fig. [21 Since there is no place for the D8 
and D8-branes to end, they are connected aiU = Uq. We also impose the boundary 
condition Aq{oo) = /i at the both ends of the D8-D8 world- volume, where is a 
constant. We will identify this constant with the chemical potential for the baryon 
number later. Solving eq. ( 13. 4p with this boundary condition we find for U ^ Uq 



This solution is singular at f/ = Uq and does not actually satisfy the original equation 
( 13. 3p unless C = 0. Therefore, we must choose C = and obtain Aq{U) = fi. 
Because of the connected configuration of the D8 and D8 branes the chiral symmetry 
U(A^/) L X U(A^/)_R on the probe D8-D8 pairs is always broken to a diagonal subgroup 
U{Nf)v in the low temperature phase. The situation is the same as in the cases 
without the gauge field background [TTl [13] . 

Instead of using the constant Uq to parametrize the solution we can also use the 
U = oo asymptotic separation L between the D8 and D8-branes in the r-direction. 
It is related to Uq by 



where t'{U) is given in eq. (13.41) with C = 0. 

Substituting eq. (13. 4p with C = into the action (13. 2p and introducing new 
variables u = U/Uq, ukk = Ukk/Uq and f{u) = 1 — u^^j^/u^ the effective action 



Ao{U) ~ Ao{Uo) + const, x C \U - f/o|i 



(3.5) 




(3.6) 
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becomes 

/■oo 

Sd8 = Ts dnu'^il 7 — — -, (3.7) 



where 

f^^M^^R^U',)^ /dV (3.8) 
9s J 

Note that this reproduce the result in ref. 



4. High temperature phase 



In the high temperature phase the geometry (12.41) is dominant. Using the ansatz 
(12.71) the induced metric on the probe D8-branes is 



^ ) ( fiU)dtl + Sijdx'dx^ 



+ 



dU' + 



u^dnl 



and the effective action of the D8-branes (12. 5p becomes 

3 



9s 



This action leads to equations of motion for r{U) and Aq{U) 

d U^f r' 

dU 




0, 



0, 



(4.1) 



(4.2) 



(4.3) 



which can be easily integrated once as before. As in the case without the gauge field 
[T3l [H] there are two types of solutions in the high temperature phase. 

One solution is similar to the one in the low temperature phase. The integration 
of eq. (14. 3p gives 



{lYf{U)(U^f{U)-Ulf{U,] 



Ao{U) = /i. 



(4.4) 



(a) 



(b) 



Figure 3: D8-D8-brane configurations in the high temperature phase. 

where Uq is an integration constant. As before we have imposed the boundary 
conditions t'{Uo) = oo and Ao(oo) = fx. A typical configuration of t{U) is shown in 
Fig. [3] (a). The chiral symmetry U(A^/) l x U(A^/)r is broken to a diagonal subgroup 
\J{Nf)v- Substituting eq. fl4.4p into eq. fl4.2l) the effective action becomes 



and Tg is given in eq. (13.81) . 

Instead of using Uo we can also use the asymptotic separation L in eq. (13. 6p to 
parametrize the solution, which is more convenient when comparing this solution to 
the other one. The relation between L and Uq is obtained from eqs. (13.61) and (14. 4p 
as 




(4.5) 



where we have rescaled the variables as m = U/Uq, ut = Ut/Uq, f{u) = 1 — Uj^/u^, 




(4.6) 



where 



F{ut) = 2 / 



oo 



du 



/(I) 



(4.7) 



Jl 




) 



For the other solution the first integration of eq. (14. 3 p gives 



r'iU) = 0, 



{2na'A',{U)f 



fj2 



(4.8) 



t/8 (f)^ + C2' 
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where C is an integration constant. t'{U) = is the trivial solution of (14.31) . A 
typical configuration is shown in Fig. [3] (b). It describes a situation that the probe 
D8 and D8-branes separately extend along the [/-direction in straight lines. The 
separation between the D8 and D8-branes is chosen to be the same as the asymptotic 
separation L in the previous solution. The chiral symmetry U(iVj)j^ x Xj^Nf)^ is 
unbroken in this case. Substituting eq. (14.81) into eq. (14. 2 p and using the rescaled 
variables as in eq. (14.51) the effective action becomes 



D8 



du 



u 



(4.9) 



where 



C2 



(4.10) 



To determine which of the two solutions is dominant we compare the values of 
the effective action. From eqs. (14.51) . (14.91) we obtain the difference as 



cU _ ell 
_ ^D8 ^D8 



du u 



u^fiu) 



u^f{u)-f{l) Vu^T^ 



du 



u 



VuF+ 



, (4.11) 



For AS* < the curved configuration (14. 4p is dominant and the chiral symmetry 
is broken, while for AS > the straight configuration (14. 8 p is dominant and the 
chiral symmetry is unbroken. Although the integrals in eqs. (14. 5p . (14. 9 P are divergent 
a.t U = oo, the difference is finite due to the same asymptotic behaviors of t{U) 
and Aq{U). We evaluate eq. (14. lip by numerical calculations. For that purpose it 
is more convenient to change an integration variable to z = u~^, which has a finite 
interval < z < 1 for 1 < m < oo. The result of the calculations is shown in Fig. |H 
The behaviors of AS* as a function of ut for various values of c are given. The special 
case c = reduces to the result in ref. [13]. In this case AS" is positive for ut larger 
than a certain value uto and negative for ut < uto- The chiral symmetry is broken 
for Ut < Uto and unbroken for ut > uto- The point ut = mto is a phase transition 
point. This phase transition is of the first order since two different configurations 
in Fig. [3] are possible at the transition point. As c increases, the transition point 
Uto decreases. When c > 0.2158, there appears a new region near = in which 
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Figure 4: AS* as a function of ut for various values of c. From the bottom to the 
top each hne represents the case for c = 0, 0.2158, 0.2252, 0.3 respectively. 



AS* > 0. When c > 0.2252, AS* is positive for all values of ut and the chiral 
symmetry is always unbroken. 

From these results we can draw a phase diagram in the c-ut space as shown in 
Fig. O The chiral symmetry is broken in the region of small c and small ut and 
unbroken outside of it. Note that we are considering here the high temperature 
phase of the confinement /deconfinement transition and only the part ut > ukk of 
this diagram is valid. 

It is more appropriate, however, to draw it in the space of the temperature T 
and the baryon number chemical potential /j,. From eq. (12.41) the temperature T is 
related to ut as 

where we have used eq. (14. 6p . 

The relation of the chemical potential /i to ut and c can be obtained as follows. 
From eq. (14.81) the large U behavior of Ao{U) has a form 

Ao{U)r^fi + ^, (4.13) 

where n and v are constants. We have chosen the same value /i for the constant 
term as in the curved solution (14. 4p . According to the AdS/CFT dictionary [1] 
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Figure 5: The phase diagram in the c-ut space. 

for a massless vector field in a six-dimensional bulk, /i is a source coupled to an 
operator of dimension four O4 on a five-dimensional boundary. The U(l) gauge 
field Aq defined on the whole of the D8-D8 world-volume contains the gauge fields 
for both of U(l)y and 11(1)^ in the fiavor symmetry (12. 6p . The part of Aq which 
is symmetric for an interchange of D8 and D8 corresponds to U(l)v', while the part 
which is antisymmetric corresponds to U(l)^ [TT1122]- Since the constant term /i is 
symmetric, it is a background value of the U(l)y gauge field coupled to the baryon 
number density O4, and fi is the baryon number chemical potential. 

Integration of eq. (14. 8 p determines Aq{U) up to a constant term (/i in eq. (14.130 ). 
We can fix this constant term as follows. We first require that Ao{U) vanishes at 
U = Ut because of the regularity. To see this we first change the coordinates from 
{U,tE) to {r,9) defined by 

f/^ = t/| + ^/^r^ (^=l{^ytE. (4.14) 

From the induced metric (14. ip with t'{U) = we see that (r, 6) are the polar 
coordinates near the point U = Ut- The point U = Ut corresponds to the origin 
r = and should be treated with care since the polar coordinates are not good 
coordinates near the origin. It is better to use the Cartesian coordinates 

y = rcos6, z = rsm6. (4-15) 

The relation between Aq and the components Ay,Az in the coordinates {y,z) is 
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obtained from AqcHe = Aydy + Azdz as 



Since we require that Ay and are regular at the origin r = 0, v4o(f/) must vanish 
a.tU = Ut- We also note that although Ao{U) is a gauge dependent quantity, it must 
vanish at U = Ut in any gauge. Only the gauge transformations which preserve the 
condition Ao{Ut) = are allowed. 

The vanishing of Aq{U) aX U = Ut fixes the constant term in this case and we 
find 

The chemical potential /i is obtained as the asymptotic value for [/ = oo 



-^{F{uT)f / duJ^ , 



/^ = ^o(oo) = — -— (FK)r / (4.18 



where we have used eq. (14.61) to eliminate Uq . This gives an expression of the chemical 
potential in terms of Ut and c. 

Using eqs. (14.121) . (I4.18P we can convert the phase diagram in Fig. Oto that in 
the /i-T space by numerical calculations. Using dimensionless variables 

T = LT. A = ^f.. (4.19) 

the phase diagram in the jl-T space is shown in Fig. El Only the part of this diagram 
for the high temperature phase of the confinement /deconfinement transition, i.e. 
T > Tc = LMkk/ i'^T^) is valid. Therefore, our result of the phase diagram looks like 
Fig. [Has we explained in Introduction. The orders of T and /i at the transition points 
can be estimated from eq. flCTD . Using = glj^NJl/ {2Mkk) and /i, T = 0(1) 
we obtain 

T = 0(L-i), iJi = 0{glj,jN,L-^M-^j,). (4.20) 

If we assume L = 0{M^]^), the transition temperature is of the order of the com- 
pactification scale Mkk, and the chemical potential is of the order gyM^cMxR, 
which is much larger than M^k since Qym^c ^ 1- 
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Figure 6: The phase diagram in the fl-T space. 

5. Discussions 

We analyzed the Sakai-Sugimoto model at finite temperature and finite baryon 
number chemical potential. The chemical potential is introduced as an asymptotic 
value of the U(l) gauge field on the probe D8-D8-branes. Using this model we 
studied the phase structure of the chiral symmetry breaking and obtained the phase 
diagram in Fig. [H This phase diagram should be compared with that expected 
in QCD [23]. Our result is especially different from the QCD expectation at low 
temperature. In QCD the chiral symmetry is expected to be unbroken even at zero 
temperature if the chemical potential is sufficiently large. In our analysis the chiral 
symmetry is always broken below since the geometry of the t-U space allows only 
the curved configuration of the probe branes. Since we used the probe approximation 
for the D8-D8-branes, the gauge field on the branes (the chemical potential) does 
not affect the bulk geometry. It is interesting to see whether back-reactions of the 
gauge field on the geometry of the t-U space change the phase structure below the 
temperature Tc. To fully understand the phase diagram we need an analysis beyond 
the probe approximation. 

In the usual field theoretical approaches to the chiral symmetry breaking one 
considers condensations of quark bilinears iptp as an order parameter. The quark 
masses are sources of these operators. In fact, in other models of the holographic 
QCD [TJ [8] the mechanism of the chiral symmetry breaking is different from that 
in the Sakai-Sugimoto model. The chiral symmetry is realized as the rotational 
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symmetry of the probe branes in the transverse space. The quark masses and the 
quark condensations can be read from asymptotic behaviors of the probe branes. In 
the Sakai-Sugimoto model quarks are always massless since the asymptotic distance 
between the D4-branes and the D8-D8-branes, which is proportional to the quark 
mass, is zero. It is not clear how to introduce quark masses in this model. It is 
interesting to clarify the relation between the mechanisms of the chiral symmetry 
breaking in the Sakai-Sugimoto model and in other models. 
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